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ABSTRACT 


The contributions of the self-phase modulation and cross 
)hase modulation, the two optical fiber non-linearities, to the 
Interference power spectrum for an intensity modulated sub-carrier 
aultiplexed, multiple optical carrier systems have been analysed. 

f 

The equations for interference power spectrum has been derived for 
i two optical carrier case and then generalized to N-optical 
carriers under some specific cases and assumptions. The signal to 
interference ratios have been calculated and plotted against the 
laser spectral width. The degradation of the above system due to 
increase in the number of optical carriers is also shown. 
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CHAPTER - 1 


INTRODUCTION 

Optical fiber haa many excellent featurea as a transmission 
medium. Among those a wide bandwidth availability is the most 
promising. Since the optical carrier frequency is in the range of 
10^^ Hz even 0.1 % utilization , of the bandwidth is sufficient to 
include ten thousand TV channels of 5 MHz bandwidth with a spacing 
of 5 MHz between them. The enormous bandwidth capacity of the 
fiber optic link is tailor made for multiple access light wave 
systems, which inc lude local area networks and metropolitan 
networks . 

1.1 SCHEMES IN MULTIPLE ACCESS SYSTEMS 

The two types of schemes which are used in multiple access 
systems are : 

a) Time division multiplexing (TDMA) 

b) Frequency division multiplexing (FDMA) 

Multiple access systems differ from wide-band links in that 
each user generally requires only a small fraction of the total 
data throughput. Uhile using the TDMA, each receiver must' receive 
all the transmitted bits and select the appropriate bits for its 
use. This unnecessarily calls for receivers with greater 
sensitivity than that^ is required- 
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The FDMA, an alternative scheme for bandwidth utilization, 
can be classified into two types : 

i) Optical frequency division multiplexing (OFDM) 

ii) Sub-carrier multiplexing (SCM) 

In the OFDM, the channel selection and other functions are carried 
out at optical frequency. This requires that the laser line width 
be very small and operated only in single mode conditions. At 
present this method of implementation is not cost effective. 


In SCM system, as shown in Fig. 1.1, data from each channel 
are used to modulate microwave subcarriers, which 
are then used to intensity modulate a single optical carrier. 
Fig. 1.2 shows optical carriers being modulated by a different set 
of subcarriers, and the signals from each set being combined in a 
power combiner and then sent in a single fiber. 

1.2 AN OVERVIEW OF OPTICAL FIBER NON-LINEARITIES IN SINGLE MODE 
FIBER 


In its simplest form, an optical fiber consists of a central 
core surrounded by a cladding layer whose refractive index is 
slightly lower than that of core (step index fiber). The two 
parameters which characterize the fiber are the normalized 
difference of refractive indices of the core (n^^) and cladding 
(n2), [1] 



( 1 . 1 ) 


and the normalized frequency V defined by 



POWER COMBINER POWER COMBINER 
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a is the core diameter 


( 1 . 2 ) 


and X is the wavelength of the transmitted wave. 

It is shown that if V < 2.405, the fiber supports only a single 
mode. The typical values of the core diameter are 25~30 ^m for a 
multimode and 2-4 pm for a single mode fiber. 

1.2.1 Attenuation 


An important design parameter is the power loss during 

transmiasion of optical signals inside the fiber. If P is the 

o 

power launched at the input of the fiber of length L, the 
transmitted power is given by 

^o (1-3) 

a = attenuation constant expressed nepers/km. 

Material absorption and Rayleigh scattering are the major factors 
contributing to fiber loss. The fiber loss at present is very 
low, about 0.12-0.15 dB/km at X = 1.55 pm. 

1.2.2 Dispersion 

In a single mode fiber, the only dispersion which is of 
importance is the group velocity dispersion. This is caused when 
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different frequencies travel with different velocities in the 
fiber. The propogation constant ft is given by 


ft 


nw 

c 


(1.4) 


The refractive index is not a constant but a function of 
frequency. This relation is given by Sellemier’s equation [1] 


nC^o) 




I 

j = l 




(1.5) 


6ij is the resonance frequency at which medium absorbs the 
optical radiation 
is the strength of resonance. 


The group velocity is the velocity with which the envelope of the 
optical field moves in the fiber and is defined by 


Vg = (^l/(d/9/di^)j (1.6) 

Substituting eqs . (1.4) and (1.5) in eq. (1.6), we see that v 

& 

is dependent on frequency. This type of dependence causes 
chromatic or material dispersion. A measure of variation of 
group-velocity with frequencies is given by the parameter ft^f 
i.e., group velocity dispersion parameter [1], 


^2 




(1.7) 




du 


where 
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1.2.3 Self Phase Modulation (SPM) and Cross Phase Modulation 
CXPM) 

As seen from the above discussion the variation of 
refractive index with frequency gives rise to chromatic 
dispersion. Another kind of non-linearity is the variation of the 
refractive index due to the incident intensity of the light. This 
effect is known as the self phase modulation (SPM), and now the 
refractive index of the fiber becomes 

Z (io, E(t)) = n(6i) + n^lECt)!^ , (1.8) 

2 

where |E(t)| is the instantaneous optical power, 

E(t) is electric field in volts/metre, 

n(«) is the linear part of refractive index and n^ is the 
non-linear refractive index coefficient. 

The intensity dependence of the refractive index leads to a 
large number of interesting non-linear effects. As evident from 
eq. (1.8), the instantaneous phase is given by 

^(t) = [n(*^) + n2|E(t)|^j L (1.9) 

L = fiber length 

X = wavelength 

= non-linear refractive index coefficient 
n(<ii) = refractive index 

The non-linear phase shift due to SPM is 

*nl. ' L f (|ECt)l2) 


( 1 . 10 ) ^ 
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The cross phase modulation (XPM) refers to the effects of the 
intensity of two copropagating fields say Ej^(t) and £ 2 ^^^ 
phase of the other. If E^(t) and E 2 Ct) are the envelopes of the 
two copropagating optical fields then the non-linear phase shift 
for the first field is given by 

LClEj^(t)i^ + 2lE2(t)}^) (1.11) 

and for the second field is given by 

^jjL(t) = K^LClE^Ct)!^ + 2iE^(t)|^) (1.12) 

where K = 2n/X. , is the wavelength of first wave, 

O 11 

= 2n/X^, X^ is the wavelength of second va.-ve. 

It can be that the XPM factor for two equally intense electric 
fields is twice as compared with that of the SPM factor. 

1.3 DETECTION SCHEME 

Ue can broadly classify the detection schemes used in 
lightwave communication into two types. 

a) Direct Detection (DD) 

b) Coherent Detection (CD) 

It is shown that the signal to noise ratio for DD system is 
lower at least by a factor of eight as compared to that of CD 
system. This advantage, however, is at the cost of receiver 
complexity and increase in the coat of a CD receiver. In a 
sub-carrier multiplexed intensify modulation system with direct 
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detection, the components required for detection are a photo-diode 
followed by a low noise amplifier and a microwave receiver. The 

microwave receiver components are easily available and 
standardized. The source line width in a DD system is not as 
stringent as in the case of CD systems. In the present work the 
intensity modulated system with direct detection is considered. 

1.4 MODULATION SCHEME 

Modulation schemes can be broadly classified into : 

1) Intensity Modulation (IM) (analog/digital) 

2) Coherent Modulation (CM) (analog/digital) 

a) Frequency Modulation (FM) (analog/digital) 

b) Phase Modulation (PM) (analog/digital) 

Intensity modulation la a simple modulation scheme to 
implement. By varying the laser bias current about a operating 
point the output power of the laser can be modulated accordingly. 
This can be detected by a direct detection scheme. 

Phase modulation and frequency modulation require that the 
laser spectral width be very narrow and the photo detector at the 
receiver has a flat response. These factors increase the coat of 
the system when we are considering multiple access system with 


coherent modulation schemes. 



1.5 OBJECTIVES AND SCOPE OF PRESENT UORK 


The effect of SPM and XPM on a single optical carier IM 
system is of little consequence, due to the inherent nature of 
square law detection, where the envelope can be recovered easily. 


Consider the IM wave with mCt) being the message and 
frequency in optical domain 

xCt) = (l+m(t))^^^ cosC^iij^t) (1.13) 

Due to the SPM, the signal would become 

x(t) = (l+m(t))^'^^ cosC^j^t + KCl+m(t))) (1.14) 

K being SPM factor. During the detection process the signal is 
squared. 

(x(t))^ = (l+m(t)) j^l/2 + 1/2 cos(2i^j + 2K(l+m(t)))j (1-15) 

^ . (l+|(t) _)^ Cltm(t ) ) ^ 2K(l+m(t)))j (1.16) 


The second term is filtered out and the message process m(t) is 
easily recovered. Now consider a two optical carrier system. The 
Individual signals are combined and sent in the optical fiber. 
The composite signal is given by 

x(t) = J 

It can be shown that after being passed through the photo-detector 
the signal has 3 components; they are 


(l+m(t)) 


1/2 


cos(oij^f ) 


]4 


(l+m2(t) ) 


1/2 


cos(<uj2f ) 


(1.17) 
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, and 

respectively (1.18) 

Not only the message processes m^(t) and m 2 Ct) are recovered but 
also the cross term is centered at | 6 i^-£ij 2 j. Without SPM and XPM 
taking into consideration, the phase factor ^(t) in the cross term 
goes to zero. It becomes evident that the SPM and XPM tend to 
spread out the power spectral density of the cross term giving 
rise to more interference. 

The purpose of this work is to analyse the amount of 
Interference caused in a multiple optical carrier, sub-carrier 
multiplexed intensity modulated system, due to self-phase 
modulation and cross-phase modulation, the two non-linearities of 
the fiber. 

1-6 PRACTICAL APPLICATION 

Consider a scheme where N customers are received by a local 
exchange, as shown Fig,. 1.4 [2]. 

The exchange provides customers’ services such as voice, 
data, and video on different sub-carrier frequencies. When the 
customers want to communicate with the exchange, each customer 
transmits the required message by modulating at a different 
microwave sub-carrier. The detector then detects the sum of all 
the optical signals. This mixing process gives rise to 
interference noise, and ultimately can limit the number of optical 
carriers and the transmission bandwidth of the message source. 




1:N COUPLER 


Figure 1.4: A Bidirectional SCM Lightwave System 
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The limitation posed by such a mixing process is analysed in this 
thesis . 

1.7 ORGANIZATION OF THESIS 

Chapter 2 discusses about the SPM and XPM factors, the 
laser phase noise, and problem formulation taking into effect SPM 
and XPM factors and laser phase noise. 

Chapter 3 deals with derivation of the cross power spectral 
density for specific cases. The signal to interference ratio 
(SIR), is calculated for various line widths of the laser source. 
The variation of SIR with number of laser sources used is 
estimated. 


Chapter 4 discusses the inference and scope for future work. 



CHAPTER - 2 


PROBLEM FORMULATION 

2.1 CALCULATION OF SPM AND XPM FACTORS 

As seen earlier the non-linear phase shift depends on the 
intensity of propagating radiation, vide eq. (1.10). If we 
consider N optical carriers travelling in the fiber with intensity 
modulated envelopes of Ej^Ct), E^Ct) E^Ct) ... Ejg(t) then the phase 
change in channel j is given by 

N 

= V L(eff) l^jE^Ct)!^ + 2 2 jE^(t)j^j (2.1) 

m =1 
m jrj 

v is called the non-linearity parameter [1] 

and w = ^ (2.2) 

ef f 

= frequency of propagating li .ght wave in radians/sec 

A = effective area in fiber from which the maximum power 
et t 

is transmitted 

2 2 
A^^^ is of the order of 10-20 ^«m in visible region and 50-80 #jm 

in 1.5 ijm wavelength region. the non-linear refractive index 

coefficient is 

n^ = 3.2 X lO"^^ cm^/U 

L = — exp (-oiL) 

eff oi ' 


(2.3) 
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a. = attenuation constant = 0.046 km ^ (0.2 dB/km). If total 

length L of the fiber is large, then 

^ l/« (2.4) 

V varies between 2-3 U ^ km ^ as the frequencies vary between 
visible and infra-red region. 

Ue define a self-phase modulation factor K 


as 

K = (2.5) 

If we assume L » l/oi = 1/0.046 km and v varies between 2-30 U ^ 
eff 

km ^ then K varies between 43.47-652.17 The cross phase 

modulation factor is twice that of SPH. 

2.2 ANALYSIS OF LASER PHASE NOISE 

The power spectrum of the laser source is not delta function 
but is considered as a lorentzian shaped power spectrum. 

The ideal laser source is represented as 

x(t) = A exp (j<ij^t) , (2.6) 

r^(x(t)) = A 6(6i - 6i^) , (2.7) 

where [ ] represents the Fourier transform operator. 

To incorporate -the spread in the power spectrum, the laser 
source is represented as 
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x(t) = A exp + ^(t)J (2.8) 

where <jJ'Ct) is usually modelled as a continuous path Brownian 
motion process with zero mean. The characteristic of Brownian 
motion stochastic process is that, the incremental processes |^(t) 
- 0(t + t)| are independent random variables with zero mean and 
with gaussian distribution [3]. 


From the Ui ener-Khint chine theorem we know that the Fourier 
transform of the autocorrelation function yields the power 
spectral density (PSD) [ ], 


PSD = (t) x(t + T)J.j (2.9) 

= I A exp - ^(t) + ^(t + ’•))jjj (2.10) 

= exp O^T) * I exp (j(<^(t+T) - ^(t))jj (2.11) 

= exp (j^-i^T) exp (-fi|T|/T^)j (2.12) 


where t is called the coherence time, and * denotes the 
c 

expectation operator. For the derivation of eq. (2.12), refer to 
Appendix A.l. 


Carrying out the above Fourier transform 


A. 2 (fiAv) 


-1 




PSD 


(2.13) 
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where Av = called the half power bandwidth of the laser 
source. The derivation of the above equation is given in Appendix 
A. 2. 


The equation (2.13) denotes a Lorentzian shaped spectra. 


Another important result which we will be making use of is 


{ 


« ■( exp (j(<^(t) + <^Ct + T)) 


} 


(2.14) 


we can write the above equation as 


{ 


* i exp (j(<^(t) - ^(t + T))+ j2^(t + T)) 


} = 


(2.15) 


Since the characteristic of Brownian process is that, the 
processes j</>(t) - ^(t+T)| are independent, we can write eq. (2.15) 
as 


{ 


* { exp (j(<^(t) - <^(t + T)) 


} *{ 


} = 


exp ( j2fl^(t + T ))i = 0 (2.16) 


It is usually assumed that the underlying process is a gaussian 
zero mean non-stat ionary random process [4]. The variance of the 
gaussian random process increases linearly with time. In such a 
situation the integral for t >> 0 goes to zero. 


i . e . 


{ 


* 4 exp ^j2(^(t + T )) 


} 


= 0 , 


(2.17) 


which is derived in Appendix A. 3. 
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2.3 SPECTRAL ANALYSIS - MIXING OF OPTICAL FIELDS 

2.3.1 Mathematical Preliminaries 


Let the real signal 


yCt) = ACt) cos (2n vt + <jfr(t)) 


(2.18) 


be a random stationary narrow-band optical process. 


Let y^(t) = A(t) exp (j<^Ct)) be the complex envelope of the 

complex signal. The subscript c denotes a complex signal. 


yj,Ct) = y^(t) exp(j2nvt) = A(t) exp ( j<ji' ( t ) ) 6xp(j2nwt), (2.19) 


which is centered around the optical carrier frequency v and is 
tightly band limited around the optical carrier frequency. The 
real signal y(t) is 


y(t) 



( 2 . 20 ) 


By the Ueiner-Khintchine theorem the power spectrum of Stochastic 
process is the Fourier transform of its autocorelation. 


The autocorrelation of the complex signal is 

#v 

where r^(T) is the complex envelope of the autocorrelation. 

The power spectrum of the complex signal is S^(t) and it is 


(t+T) y* (t)) 


} 


exp ( j 2nvT ) = r (t ) exp ( j 2nvT ) , 
c 


( 2 . 21 ) 


given by 
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Sc(t) = <^2.22) 

To find the power spectrum o£ the real signal we have to find its 
autocorrelation function and then take the Fourier transform of 
the autocorrelation function, 

(t)j] (2.23) 

Since, it is known that. 

Re { }Re { = 1/2 Re { AB^} + 1/2 Re { AB* > (2.24) 


r(t, t+T) ='&r|Re |^y^Ct+T)j. Re-j^: 


we have 


E{ycCt+T) y*(t)| 


r(t, t+T) = 1/2 Re *^y^(t+T) y^,(t)J + 1/2 Re « -(ycCt+T) y^(t) 




} 


(2.25) 


r(t, t+T) = 1/2 Re {r^(T)) + 1/2 Re*£{<y^(t+T )"'y^(t) exp (j2nv(t+T):0 


(2.26) 


where <.> indicates time average. 


Since we are interested in the average power spectral density we 
have to do a combination of statistical and time averaging. It is 
shown in [4] that the second term of eq. (2.26) vanishes. 
Therefore, the autocorrelation of the real signal is half the 
real part of the complex signal. 


r (T) = 1/2 Re { T^(t) } (2.27) 

c c 

2.3.2 Problem Definition 

The general equation for the electric field of an aonfUtuclt 


modulated signal is 
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(2.28) 

1/2 

where (p^Ct)) is the modulated envelope and is an optical 

carrier frequency. If there are two optical carriers the 

the combined equation of the electric field is 

E^(t) = Cpj^(t))^'^^ exp (j^^j^t) + (p 2 (t))^^^ exp (2.29) 

1/2 

where (P 2 C't)) is the modulated envelope and is the second 

optical carrier frequency. 

If the fiber optic link does not introduce any distortion, 
then the signal incident on the photo detector would be 

[ E(t) = [ Re { E^^(t) } + Re { (2.30) 

At the output of the photo detector we would have 

i(t) = a: (E(t))^ (*) h^(t) , (2.31) 

where x is the responsivity of photo detector, amper es/watts 
i(t) is the output current 

h^(t) is the impulse response of the photo detector 
(*) denotes the convolution operation. 

For the purpose of analysis, we assume that the photo detector be 
an ideal low pass filter with a cut off frequency which is much 
lower than the optical frequency. The resulting photo current 
after low pass filtering of eq. (2.30) yields [5] 
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iCt) = P^(t) + p^Ct) + 2 Re { (p^Ct) P 2 Ct))^/^ exp C J ^ ^ 

= i^Ct) + i2(t) + ij^(t) (2.32) 

The first two terms ij^(t) and i^Ct) give us the modulated 
envelope. The third term ijj.C't) however, is the result of mixing 
process and represents an undesired interference component. If its 
power spectrum lies within the bandwidth of the power spectrum of 
the envelopes ij^(t) and i^Ct). 

To find the power spectral density of the resulting photo 
current i(t), we must consider the Fourier transform of the 
autocorrelation, i.e. 

* [ i(t+T) i*(t)j j = S(f) (2.33) 

Uith the assumption that Ej^^(t) and zero mean 

independent processes, it is shown that S(f) can be written as [5] 

SCf) = * S 3 CO * {r^TT^j } C2.34) 

where, 

Sj^(f) is the power spectrum of Pj^(t) 

S 2 (f) is the power spectrum of p 2 (t) 

S^^CO ia { S Ci^(t)) * » (i2Ct)) > » 

ij(t) and l 2 Ct) contain sinusoidal modulation terms. The Fourier 
transform shown above yields the impulsive cross power spectrum. 


II I 2 *5 (f) , 
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where and 1 2 are two constants representing the average DC 
values of i^Ct) and i 2 Ct). Note that rg 2 denotes the time 

averaged autocorrelation of the signals. Since Ej^(t) and ^ 2 ^^^ 
are independent, then ^£^^^£2 s^P^^rates into F^^ ^E2 ' they 

are time separable [5]. Therefore, 

■17^2' ''eI ''e 2 = * { ^itt) ] * { EjCt^) EjCt) I 

where Ej^(t) = Re { E^^(t) } 

= Re { Ej^^(t) exp (j2n i^^t) } 
and ^2^^^ ” ^ ^ 

= Re { E 2 ^(t) exp (j2n V 2 t) } . 

Note that the over bar is used to denote time averaging in the 
above equations. 

a) The Effect of Adding Laser Phase Noise 

As shown previously, the laser phase noise would add a extra 
factor exp (j<j^j^(t)) to the signals, i.e. consider 
E- (t) are two independent optical fields with a modulated 

2 * Cl 

envelope 

Eic(t) = (Pj^(t))^^^ exp exp (2.35) 

and 

E 2 g(t) = CP 2 (t))^^^ exp (2.36) 

Now the complex envelopes Ej^^(t) and E 2 j,(t) are 
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6XP (j4!>j^(t) 

(2 

.37) 



exp (j^!/ 2 (t) 

(2 

.38) 


After squaring of the combined fields at the photo detector the 
cross term can be represented as 

i^(t) = 2 Re CJ2n-v^tD Xexp (.-j2nv^t:)^ (2.39) 

Note from eq. (2.34) that the terms contributing to interference 
power spectrum is nothing but the Fourier transform of the time 
average autocorrelation of the cross term. 

Let and F^^ represent the autocorrelation of Ej^(t) and 

E 2 (t), then 

■■eI '■e 2 = 1/2 E^^CtJ E^^Ct)] j 

+ 1/2 Re ^ E^^(t+T) E 2 *(t+T) E^^(t) E 2 *(t)J J (2.40) 

To calculate the interference power spectrum, we should now 
calculate the average power spectrum of the first term, i.e. 

r^l^l /2 Re ^ [pj^(t+T) p^(t) P2(t) P2(t+r)j 

X exp (j(<^^(t+T) - </ 2 j^(t)) - K4’ ~ (j2nAvT)jj 

(2.41) 

where, 

difference frequency. 
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<jfr^Ct+T) - is the incremental phase noise process due to 

laser source I, 

<(^ 2 (^+ 1 ’) - is the incremental phse noise process due to 

laser source II, and 

p^(t) and P 2 *'^^ envelopes of the two signals. 

(b) To Account for SPM and XPM Effects on Signal 


Now the signal leaving the optical fiber has to be further 


refined 

to consider to 

eff ect 

of SPM and 

XPM. As before we 

consider 

complex signals 

Ei,Ct) 

and E2^(t). 

The incident 


intensity on the photo detector would be 


[ E(t) = (Re { E^^(t) + Re E 2 ^(t) 

where ^*^4 E 2 ^(t) are given by 


E^^Ct) = (p^Ct))^/^ exp(j6>^t + jK^p^(t) + j2Kj^P2(t) + j>j^(t)) 

(2.42) 

E2cCt) = (p2Ct))^'^^ exp(j6i2t + jK2P2(t) + j2K2p^(t) + j<}^2Ct)) 

(2.43) 

n^ <<>- L .. 

and, = SPM factor 1 = = ■ ■ - (2.44) 


2K^ = XPM factor 1 


(2.45) 


K 2 = SPM factor 2 




“2 "2 ^e£f 


C A 


eff 


2K, 


(2.46) 


XPM factor 2 



pj^(t) and p 2 Ct) contain the message processes which may represent 
any type of modulation scheme, and are laser phase 

noise of sources I and II respectively. 

The above equation can be written conveniently as 

EicCt) = 6Xp(jv»^(t)) exp (jojj^t) = Ej^^(t) exp 

(2.47) 

E 2 c('t) = (P 2 ('t))^^^ expCjv^Ct)) 6xp (j^^ 2 ‘‘^^ “ ^ 20 ^"*^^ 

(2.48) 

where 

V^(t) = K^p^(t) + 2K^P2(t) + ^^(t) 

V'2Ct) = K2P2(t) + 2K^p^(t) + 

Since Pj^Ct) and P 2 C't) contain the message processes, <^j^(t) and 
^-(t) which are random phase fluctuation processes, we observe 

4b 

that the signals Ej^^(t) expCji.>j^t) and ^2c^^^ 

random narrow band optical processes centered around the two 
optical frequencies and respectively. 

Uith E. (t) and E- (t) are given by 
1 C A c 

El^jCt) = (p^Ct))^/^ exp (jKj^Pj^Ct) + j2K^P2(t) + j<^j^Ct)) (2.49) 

E2c(t) = (P2Ct))^''^ exp(jK2P2Ct) + j2K2P^(t) + 3-5^2^^^^ (2.50) 

as shown earlier, the resulting photo current is given by 
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where ;t is usually taken as 1 amp/watt. The resulting current 
could be written as 

i(t) = Pj^(t) + P 2 C‘t) + 2 Re { Ej^^(t) E^^Ct) {expCjZnAvt) } 

(2.52) 

The total power spectrum o£ the photo current would be 

S(f) = [* [ i(t + T) i (t) jj . (2.53) 

As seen from earlier discussions, the significant contribution to 
the cross power spectral density S^(f) is the power spectrum of 
the cross term. 

S^(f) = ^x[ * { E2*(t+T) j Re[E^*(t) E2(t)j J j 

(2.54) 

The main thrust of this thesis is to derive and evaluate eq. 
(2.53), the cross power spectrum under some specific cases. 

After evaluting eq. (2.53), we can find the signal to 
interference ratio, SIR, as 

_ (Total message power) 

(Interference power within in message bandwidth) 

J S (f) df 

SIR = = (2.55) 

f S (f) dt 
J X 

where the integration of S^(f) is carried out over message 
bandwidth and, Sj^(f) is the received power spectrum of the message 


PjCt). 



CHAPTER - 3 


ANALYSIS AND DISCUSSION 

3.1 THE BASIC STRUCTURE 

The block diagram in Fig. is a simple two optical carrier 
system, with direct detection scheme, and the envelopes being 
sub-carrier multiplexed. Such a scheme is analysed, and the 
signal to Interference ratio (SIR) parameter which is a measure of 
the effect of SPM and XPM is evaluated for 2 optical carrier 
system; the resultis generalized to systems with more than 2 
optical carriers. 

3.2 ASSUMPTIONS 

In subsequent analysis the following assumptions were made 

a) The only type of non-linear e^ectsconsidered are the SPM and 
XPM. 

b) The type of fiber la single mode. 

c) The photo diode acta as an ideal low pass filter with a cut 
off frequency in the microwave region. 

d) The effect of noise introduced at any stage of the system 
(Fig.2T) i3 considered negligible. 




POWER COMBINER 


Figure 3.1: A Simple Two Carrier Light Wave System 


8 
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3 . 3 CASE ONE 

At first, we will consider the case of two sine waves at 
microwave frequencies f^ and f^ intensity modulating two optical 
carriers and The optical carrier being combined in a power 

combiner, launched in the optical fiber. At the receiver the 
combined signal la passed through a photo diode and a low noise 
amplifier, the two microwave frequencies f^^ and f^ at'e recovered 
using microwave receivers whose bandwidths are centered around f^ 
and f , - Theoretical evaluation of SIR of such a case is presented 

lb 

below, wherein the effect of laser phase noise is considered. 

Let E^ (t) and E- ft) be the two complex electric fields 
Ic^ 2c' 

with w and v being the optical carrier frequencies. 

1 u 

The combined electric field E(t) 

E(t) = Re { E^^Ct) ) + Re { E 2 ^(t) } (3.1) 

where 

Ej^c^^) = (Pj^(t))^^^ exp Oj^t + jKj^p^(t) + 2 jKj^P 2 Ct) + j<^j^Ct)) 

(3.2) 

E 2 g(t) = (P 2 (t))^'^^ exp 02 ''^ 2 jK 2 pj(t) + 

(3.3) 

and, 

V -iK- ^ 

1 *“ C A ,, oi 

ef f 


SPM factor 1 


(3.4) 
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2 2 

^2 “ CT factor 2 ( 3 . 5 ) 

eff 

Refer to Section 2.1 for the definitions of and K^- 

Furthermore , 


P^ct) 

= P^Cl + 

sin(io^t) 

with [Aj^ sin(6ij^(t ) ) 1 < 1 

( 3 . 6 ) 

P2(t) 

= P2(l + 

A2 sin(6i2t) 

with IA2 sin(6>2(t))| < 1 

( 3 . 7 ) 

and, 

and p2 are 

average laser 

source powers. 



and k^ are constants and control the depth of modulation. 
For the sake of simplicity we assume 

= ^2 Pi " ^2 

The incident intensity on the photo detector would be 


[ E(t)]^ = (Re { E^^(t) } + Re { E2^(t) })^ 


( 3 . 8 ) 



+ — ^ + Re { E^^(t) £2^ (t) } 

+ Re { E^^(t) E2^(t) } 

+ Re { E2^(t) E2^(t) > ( 3 . 9 ) 

+ Re { E^^(t) E2^(t) } 


It can be easily shown that terms like Re { Ej^^(t) Ej^^(t) }, 
Re { E2^(t) E2^(t) ), Re { E^^(t) E^^(t) } have neglieible 
frequency contents in the region of operation of the photo diode. 



31 


The above mentioned terms are discussed and analysed in Sec 3.2 
under a more general case. 

The term Re { ^ centered around - ^ 2 !’ 

this would lie within the cut off frequency of the photo diode. 

On the assumption that if and are closely spaced 

optical carriers such that - i>^\ ^ 0 then S K 2 , we can 

* 

write the term Re after low pass filtering as 

> “ '‘*{22 exp Jjnwj^Ct) + jlio^Ct) 

n 1 

+ (3.10) 

The above result is derived in Appendix B.l. 

Here the coefficients C , are determined by Bessel 

nl ' 

coefficients, J (/3. ) and J-(^o)» where ft. and ft depend on SPM and 
n 1 1 z 1 z 

XPM factors, average power of the laser sources, and the factor 
controlling the depth of modulation. Further information is given 
in Appendix B.l. 

Ue know that the Fourier transform of the autocorrelation of 
eq. (3.10) would yield the cross power spectrum. 

Re (A^) Re (B^) = 1/2 Re (AB^) + 1/2 Re (AB*) (3.11) 

Taking the autocorrelation of eq. (3.10) we get 
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= 1/2 Es 2 2 2 '^nl '=n'l- “P + 

n 1 n’l’ 

+ (1 - l’)jij2t + - j<^^(t+T) + ” 



(n-n')j/o^t + 
)} ] 


1/2 Re [*{nn ‘^ni ^n’l’ + Cn+n’)j£i^t + 


n 1 n’l 


+ (1 + ■*■ + j4^l(t+T) - j4fr2Ct+T) - 


12 ) 


Since we are conaiderlng a determiniatic aignal, the expectation 
operator ia only for phase noise processes <^j^(t) and with 

the characteristic that, (refer to Section 2.2) 


X 


X 


{ 

{ 


exp (jl<^(t) -^(t+T)|) 
exp ( j|<^Ct) + ^(t+T)| ) 


exp C “ I T I nAu ) 


0 


(3.13) 


(3.14) 


where 1 /t = Av 
' c 


Aw represents the half power laser spectral 
width. 


Using the result of eq. (3.13) and (3.14) in eq. (3.12) and using 
the notion of independence between the phase processes ^j^(t) and 
02(t) eq. (3.12) can be reduced to 
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[*{e, 


^(t, t+T) E2^(t,t+T) 


}] ■•[■{ 




1/2 Re Hnn ^nl ^n’l’ CjndJ^T + (n-n’)j<.i^t + 1 


n 1 n’l 




+ (1 - j^^Ct) - j<!^j^(t+T) + j<ji2Ct+T) - ] C3.15) 

Since we ane inter'ested. in, the avenage powen apectcuni we 
perform the time average of eq. (3.15) before computing the 
Fourier transform. It is shown in Sec. 3.2 for a more general 
case that, the time averaging converts it into a double summation. 
Now eq. (3.15) becomes 

“ 1/2 Re I 2 J exp (jn^j^T + (“1'’^! 2n (AiAj^+Av^ ) 


(3.16) 


Ue can easily compute the Fourier transform of eq. (3.16). This 
derivation is given in Appendix B.2. 


S^(f) 


22 ""nl C 


nl { 


_2 

(nAv> ) 


- ^nl) 2 


2(f - f 


) 2 ] } 

-] 


C3.17) 


Where f^. = nf. + If- 
ni 12 


and, f 


= i^^/Zn and f 2 


<^2^l2n 


are the two sub-carrier 

^ j;- 


frequencies . 
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3.3.1 Expression for evaluation of SIR 

The average power of the received single tone signal is 
2 

l/4(Aj^) - Ue assume that the microwave receiver bandwidth is bHz. 
Then the expression for SIR for the microwave receiver receiving 
the single tone signal at frequency f^ is given by 


1/4 (A^)‘ 


(at £,) 


(3.18) 


^ J , S^Cf) dt 


fl-b/2 


Since the Lorentzian shaped spectrum has a smooth variation, we 
can approximate it by a constant over a small region of interest. 
Eq. (3.18) can now be written as 


1/4 (A^)‘ 


SIR = 

(at f^) S^(f^) X b Hz (B.U. of R^) 


(3.19) 


Example : 


Let 


fj^ = 2 GHz and f^ = 2.01 GHz 


Pi = P 2 = 1 mU, the average laser power 


X ^ X^ = 1.55 /Jm 


Then 


where 


= 56.398 U 


A = 50 ^im^, « = 0.046 km"^ (0.2 dB/km) 

ef f f- » 


Ai = A^ 


0.75 (75 % modulation) 
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Bandwidth of receiver S 5 KHz (centered at f^) 

= 50 MHz , The total half power bandwidth of 

laser source I and laser source II added. 


The coefficients, i were calculated and a computer 
program was written to evaluate the following equation. 


SIR(dB) = 10 log 
(at f^) 


1/4 (A‘) 


L X b Hz (B.U. of R^) J 


(3.20) 


It was found that 

SIR (dB) = 44. la 

at 

and 

SIR (dB) = 44.18 

at f^ 

3 . 4 CASE TUO 

In this section we consider the case of two data streams 
m (t) and m,(t) both being BPSK modulated with sub-carrier 
frequencies of and €2 

mj^(t) = Aj^ 3in(2fifj^t + «j^(t)) (3.21) 

m2(t) = A^ sin(2n£2t + (3.22) 


where , 
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«lCt) 

= 2nf^ 2 S(t - nT) 

j 


«2Ct) 

= 2nf^ .J act - nT) 

k 


where , 



'd = 

phase modulation index 


^k’ 

= i 1, binary data streams 


g(t) is 

unit rectangular pulse 


T is the bit period 


Now the modulated envelopes of the optical carrier 

frequencies 

and are 



(p^Ct)) 

= (PjCl + m^(t))^/2 

(3.23) 

(P2(t)) 

= (P2^^ 

(3.24) 

Pj^ and p^ are 

average laser powers. 


The two 

electric fields Ej^^(t) and E 2 ^(t) are 

combined in a 

power combiner and sent in the fiber. The incident 

intensity on 

the photo detector would be 


[E(t)]^ 

= [Re { E^^(t) } + Re { E 2 ^Ct) }]^ 

(3.25) 


where 


exp + jKj^Pj^(t) + 2 jKj^p 2 (t) + j^j^(t)) 


(3.26) 



:57 


1/2 

^ CP 2 Ct)) exp + jK 2 P 2 Ct) + ajK^p^Ct) + 

(3.27) 

In the above equations, the symbols have their usual meanings as 
explained earlier. 

After squaring operation we get 

[E(t)]^ = (p^(t))/2 + (p2Ct))/2 + Re { Ej^^(t) E 2 *(t) } 

+ Re { E^^(t) E^^Ct) } 

+ Re { E^^Ct) E 2 ^(t) } (3.28) 

+ Re { E^^(t) E 2 ^(t) } 

Consider the term Re { E^^(t) E^^(t) ). By substituting for E^^(t) 
we get 

= Re { (pj^(t)) exp (2jVj^t + 2jK^p^(t) + 4jKj^P2(t) + 2j^^(t)) } 

(3.29) 

= Re { p^(l + Aj^sin(i^j^t + aij^(t))exp (2ji>j^t + 2j<^j(t) + 4 ^^ 

+ sin(w^t + oij^(t)) + sin(i.>2'‘= ^ (3.30) 

= Re|p^(l+A^sin(6>^t+«^(t))r2 2 C/32 ^ ®*P (2>^t+2 (t) + 

+ n(f^j^t + oij^(t)) + lCi^2'*^ "2^“*^^^] } (3.31) 


By substituting. 
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sinCwj,t + oi^Ct)) "" i [ + «j^(t)]- 6xp [-jiOj^t - jot^Ct)]j 


equation (3.31) can be written as 


= Re {l 2 '"nl exp (2 jw^t + 2 j<^^(t) + + jnCi^j^t + not^Ct)) 


+ jCl^^j^t + lat2Ct)) J 


(3.32) 


where 


JP A 




(3.33) 


(refer to Appendix B.l for further details). 

For |n| > 2 and |lj >2 the coefficients C go to zero. 


Hence it is evident from eq. (3.32) that the contribution of 
frequency to the region of operation of the photo diode is 
negligible and therefore can be neglected. Similarly it can be 
shown that Re { Ej^^(t) E 2 ^(t) ) and Re { ^2c^^^ ^ have 
also negligible frequency content within the cut off frequency of 
photo diode. 


Now consider the term 

Re { Ej^^(t) E2 c^^^ ^ substituting for E^^(t) and E 2 ^(t) 


we get 
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Re 

= Re |CPi(t) FcjUj^ - 

+ jK^p^(t) + ZjK^p^Ct) 


> 2 ) t + j<^j^(t) - 


jK2P2Ct) - 2jK 


iViCtj] I 


(3.34) 

where 


Pj^Ct) = (PjjCl + Aj^sinCWj^t + «j^Ct)))) and |Aj^ sinCw^^t + a^Ct))| < 1 
PzCt) = (p^(l + A23in(<.j2^ ^ i^2 ^ ^ 

(3.35) 

The aaaumption that and are cloaely spaced, and so “ 

0 implies Kj^ « K 2 with this assumption. Uhen eq. (3.35) is 
substituted in eq. (3.34) we get 

Re jE^^Ct) E^^Ct)} 

Re-|^[[p^(l + sin(iJj^t + ciij^( t ) ) ) ] [p^( 1 + A 2 sin(6J2t + a2(t)))]j 

exp(- jKj^p^Aj^ sin(<^j^t + «^(t)) + jKj^p^A 2 sin(6ij^t + « 2 (t)) 

exp (j<^j^(t) - (3.36) 

By taking the approximation that 

(1 + x)^^^ = 1 + x/2 ; 1*1 < 1 

The envelope of eq. (3.36) can be written as 



‘♦u 


(PqC^ + sinCoi^t + ci^(t)))) (p^(l + sinCii^t + 

+ sin( 6 )^t + «j^(t)) sinCw^t + ^^(t)) (3.37) 

The complex phase part of eq. (3.36) can be expressed in terms of 
Bessel's function [ 6 ] that is 

exp ■*■ “j^Ct)) + jKj^p^A 2 sin( 6 J 2 t + 

(WjCt) - Wj'*) } C3.38) 

=n^ J^(/ 32 )C- 1 )” exp (jn [ <^^t + a^(t) ] + jl [ 6 i^t + ^^(t)] 

+ j<Ai(t) - (3.39) 

where 

= K.p A. and ft- - K^p A_ 

' 1 1*^0 1 ‘ 2 1*^0 2 

On substituting the envelope part approximation, that is, 

eq. (3.37) and the phase part expression in eq. (3.39) and 
carrying out the analysis similar to that shown in Appendix B.l, 
we can write eq. (3.36) as 

= Re ^ ®*P Cjn + oi^(t)] + jl [ 6 i 2 t + ot 2 (t)] 

+ j^^Ct) - (3-40) 


The resulting photo current can now be expressed as 



41 


ict) = 


p-(t) p^ct) 

— + Re 


{u 


^nl ■*■ 


+ jl [i^^t + oi^Ct)] + (cross term) ( 3 . 41 ) 

To calculate the cross power spectral density we have to calculate 
the cross power spectrum of eq. ( 3 . 12 ). The autocorrelation of 
the cross term is 


1/2 Re Hnn ^nl ‘^n’l’ 


+ (n-n’)j«^t + lj^ 2 ^ 


n 1 n’l’ 


+ (1 - l’)>2'*^ * jnoij^(t+T) - jn’aj^(t) + jla2(t+T) - jl’oi2(t) ^ 
j<^j^(t) - j^j^Ct+T) + j^ 2 ('t+T) - j 

+ 1/2 Re “ 1 ^ 2 2 2 2 ^nl ^n’l’ + (n+n’)j<.ijt + lj<.i. 


n 1 n’l’ 


+ (1 + l’)jii>2'*^ J*“*j^(t+T) + 

j^j^(t) + j<^j^(t+T) - 


jn’«^(t) + jl« 2 (t+T) + jl’oi^ft) + 

- 4 1 ( 3 . 42 ) 


Ue define 


y ^^ (t , t+T ) = exp (jn«j^(t+T) - jn’aj^(t)) (3.43) 

and 

yj^j^(t,t+T) = exp (jloi 2 ('t+T) - jl’« 2 C'‘^)) (3.44) 

Since the two message sources otj^(t) and ot 2 C't) are independent, 

y^^(t,t+T), y^j^(t,t+T) are also independent stochastic processes, 
(t) and <^-(t) ahe also two independent non-stat ionary zero mean 

•4 

It is to be noted that ^ 1 ^^'*' » ^ ' 


stochastic processes. 
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, t+T ) , S’!! independent stochastic processes. 
By invoking the property of independence we can write eq . (3.42) 
as 


1/2 Re 2122 '"ni Cl exp (jn<.ij^T + Cn-n’)jf.ijt + Ijio, 


n 1 n’l' 


+ (1 - I’)j6i2t) « * |^yiiCt,t+r )j. 

■j^expC ) - >;^^Ct+T ) )j. X -j^expC j<^2C^^ (first term) 


+ 1/2 Re 


[uu 

n 1 n’l’ 


Cni C^.j^ exp (jn6>j^T + (n-n’)j4i>^t + 


(1 - l')j< 02 t) « ^exp ( jnaj^(t+T)+ jn’«^(t))| 
*-|^exp (jloi 2 (t+T) + jl’oi 2 Ct))| 


*|^6xp ( j<i6j^(t) + >^^(t+T ))j- S ■|^exp(- 5 (second term) 


(3.45) 


By the aplieation of the result of eq. (2.14) the second term of 
eq. (3.45) goes to zero and only the first term remains. 


Since we are interested in the average power spectral 
density we have to compute the time average autocorrelation 
function, time average is over variable t. 


Since yuC"*- » ^ 

independent stochastic processes, the time average of the 
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first term o£ eq. (3.45) is the product of the time averages of 
each of them, and they are time separable [5]. That is, 


^ ^E1 ^E2 ^ 


< ’^El > < ^E2 > 


(3.46) 


where < > denotes the time averaging operation. 


After doing the time average 


of the first term of 


eq. (3.45), we get 






n 1 n'l' 


< S |exp (j<^i(t) - j^i(t+T) + jn6i^T + (n-n*)>^t) | > 

< * |exp(3^2^t+T) - 3 ^ 2 ^^^ j^“2^ ^ (l-l')3-2t) | > ] 


(3.47) 


Ue know that 


|exp (3<j^i(t+T) -3<^^(t))| = exp (-"ItI/t^^) 


(3.48) 


and * (axp = “I’ 

These results are derived in Appendix A.l. 

By applying tha raaulta o£ aq. d-lS) a-"! ‘I- 

follows that 

< X J exp (3(<i«'i(t+T) - ^iCt)) - 3(n-n’)^lt) > 

< exp(-n|T|/T^^) J exp(3(n-ti’) ^^t) dt C3.50. 



Since we know that 


1 

T 

+T^/2 

J expCjCn-n’) t) dt = 

- 1 

0 ; for n n’ 



-T^/2 

1 ; for n = n’ 


i 

+T 2/2 



1 

T 

J expC j(l-l’ ) tj^t) dt = 

0 ; for 1 ?! 1 ’ 



' --^ 2/2 

1 ; for 1 = 1’ 


where T^ 

= and T, = 

2 <^2 



The eq. (3. 

47) exists only for n= n 

and 1 = 1’ 

and by 

substituting 

the same eq. (3.47) reduces 

to double summation 

= 1/2 Re 1 J J C*3 v-jCt) VjCt) e*p( 

-ItI/t^i) exp(-|T| 

!'"c2 ^ 


n 1 




exp (jn 6 >^T + jl 6 > 2 T) ^ 


(3.51) 

It la shown 

in [7] that the time average 

autocorrelation 

of the 

following are functions of t only. 



< X 

l^exp Cjn(a^Ct+T) - a^(t))j- > 

= V'j^C'r) 

(3.52) 

< X 

l^exp ( jl(oi2(^''’'’^ ) “ ^ 

= ¥^ 2 ( 7 ) 

(3.53) 

The fourier 

transform of the eq. (3.51) 

is 



S^Cf) = 1/2 Re fn |C„ll^ S^(f) (*) S3Cf) (*) 

n 1 





45 


= * ^lor2^ • 

n 1 

■^22l^nll ^2^*^ ^lorl^’^ ^lor2^ • 

(3.55) 

Under- the aasumption that the half power bandwidtha of the laser 
source 1 and laser source 2 are much greater than the message 
bandwidth, the equation of S^C£) reduces to 


s^Cf) = V2 n l<:„il'[5lc [ — r 4a ■ - f 0 -.2 


+ 2Y£ I £ N 2 1 1 (3,56) 


Since a 

microwave 

receiver is used whose 

bandwidth 

is 

equal 

to 

message 

bandwidth 

at the 0/P of microwave 

receiver , 

the 

SIR, 

say 


for message one would be 

1/4 (m (t))^ 

SIR = i 

message * bandwidth of message one 


and for message two would be 


SIR 

message 

two 


1/4 Cm2Ct))^ 

S (f,) X bandwidth of message two 
X, 4* 


Thia result can easily be genralized to more than two 
carriers, consider 3 optical carriers, each of them have messages 
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p^(t), P 2 ^'^'> su-b-carrier centered at £^, and f^, 

then the resulting field incident on photo detector would be 


2 Pi(t) ^3^^^ 

[ECt)] = — 2 2 — ^ — 2 — ^ cross term (3.57) 

Assuming that the optical frequencies are closely spaced, i.e., 

= V- = V . The cross terms that would contribute to the 
1 ^ 

cross power spectrum will be 

each of them can be approximated by a bessel function, as shown in 
Appendix B.l. Theresulting cross power spectrum can be written as 


S (£) 
x' 




^23x<'> 




where 


^12x*'^^ ^ ^23x^^^ ^ ^31x^^^ are the cross power spectra 
the factors Re { E^^(t) E^^Ct) } , Re { E^^Ct) E^^Ct) }, 
{E^^(t) E*^(t) >reapectively, where ’ ^23x^^^ 

are similar to the power spectrum expressed in eq. (3.56) 


due to 
and Re 




A graph is drawn with the laser spectral width variation as 
the abscissa and the SIR variation as the ordinate. For a 2 
carrier (Fig. 3.2) and 3 carrier (Fig. 3.3) case the following 
were the values of different parameters 


^ = 1-.55 iJm 

= 53.45 U"^ 

A = A- = 0.75 (75% modulation) 

X iS 
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Average laser power o£ all lasers is taken to be 1 mU. The half 
power laser width Av is varied between 100 KHz - 10 MHz, the 
sub-carier frequencies were centered at f = 2.0 GHz, f = 2.01 
GHz and f^ ~ 2.1 GHz with a bandwidth increasing from 5 KHz to 
25 KHz in steps of 5 KHz. 

As seen from the above graph, as Av increases SIR increases. 
This can be explained by the fact that the average laser power 
being constant, increasing Av decreases the integral, 

f+B/2 

J S^Cf) d£ 

f-B/2 

But as Av is increased more, the strength of interference power 
contributed by other terms increases and, therefore SIR variation 
is not significant. 

Figure 3.4 shows the variation of the number of laser 
sources used in the system with SIR As is evident from the figure 
the SIR decreases with the number of laser sources increases. The 
spacing between sub-carriers was 50 MHz with the first sub-carrier 
at 2 GHz. The modulation index was 0.75 (75% modualtion). The 
SPM factor was equal to 53.95 U ^ and laser width was kept at 50 


MHz. 



CHAPTER - 4- 


RESULTS AND CONCLUSIONS 

It was noted that the interference noise limits the maximum 
number of optical carriers in an intensity modulated sub-carrier 
multiplexed system. The amount of interference power increases 
when the bandwidth of the message modulating the sub-carriers were 
increased. Hence the number of optical carrier used reduces 
considerably. 

The variations of the SPM and XPM factors for a fixed laser 
power did not considerably affect the SIR. There was a variation 
of about + 0.5 dB as the SPM and XPM factors were increased. It 
can also be seen that the increase in the laser power, would 
increase the coefficients, having significant contribution to 
the cross power spectral density, which would be detrimental to 
SIR. It was also noted that if the spacing between the carriers 
were reduced, the SIR did not deviate from ± 1.0 dB limits. The 
present thesis has considered only the cases when the message 
bandwidth is much smaller than the half power laser spectral 
width. 


The derivations and the consequent results were obtained 
under the assumptions that the optical carriers were closely 
spaced. This led to certain simplifications in derivations. An 
interesting problem would be to consider optical carriers sparsely 
placed. It would also be significant to consider the effect of. 
the message bandwidth in the derivation of the PSD. 


11543 
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APPENDIX - A.1 


To prove 


t I 6XpCj(^(t4T) - 
Ue state a standard result [10] 


>1 


+00 

J expC-a^b^) cos (bx) dx 
-00 


exp (-jr I /tc) 



(A. 1.1) 


(A. 1.2) 


Since <^(t+T) - ^(t) has a gaussian distributed random variable 
with zero mean and whose variance is 2n|T|/Tc . 


Ue know that if x is a random variable with a density 

function p(x) then the expectation of g(x), which is a function 

of X, is given by [8] 

+<» 

* ^ &Cx) J = J e(x) p(x) dx 
-ao 

Using above equation we get 

» ^ exp(j(^(t+T) - 0(t)) 



yjjs, /-j2,^2. +00 . ,.2,»2. 

exp Cj^fr) 6xp(0 /2& ) cos ^ exp (-^ /2& ) 

f J.f - J 


-00 


iTSicf 


*00 


■/Sno' 


i J 


2 2 

+ sin ^ expC-^ /2<y ) 




dp 



53 


2 2 

Since sin 4 > is sn odd function and (exp (-^ /Zo- )/-/2my) is 
function the second integral becomes zero and we have to 
the first term only. 


(cos 4 ^) exp (-^^/Zo-^) 


-00 




r cos(^) exp (■ 


where a 


Zo- 


Z ’ 


implies o- = 


(2a2)l/2 


^ T "/^ a X "fit 

the integral = exp 

TJn a 




= exp (- ^ /Z) 



an even 
evaluate 

2 2 ... 
a X ) 



APPENDIX - A.e 


To ahov that 


{ 


-j exp ( + Jw^T) exp (-ii|Tj/Tc) 


} 


C2(nAv)) 


-1 




Ue know that the following are tranaform paira [9] 

exp ^ 

l + (2nf)^ 

and the following propertiea are alao known [9] 


if 

e(T) ^ 

G(f) 


then 

g(aT) 

^ G (£M) 

- property 1 

and also 

j2nf 

eCt) e-* o 

^ ^ G(f - f^) 

- property 2 

where 

a = n/Tc 




Applying the tranaform pair and uaing the above propertiea 
exp (-n|T I /Tc) j 

2CnAv)"^ 



and from property 2 ve have 


exp ( + j<^^T) exp (-n|T|/Tc) 


A (2(nAv) 


-1 


r 1 + 


L 

[ nAv> J 


where = 2fi£ 

o o 

and ^ c " * half source laser spectral width. 
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APPENDIX - A.3 


To prove that 


t -I exp I j (2^^( t+T ) ) I V = 0 , for sufficiently large t. 


^Ct) is usually modelled as a continuous path Brownian motion 
process with zero mean [11]. The random variable for the noise 
process ^(t) is assumed to have a gaussian distribution [7]. The 
variance of the gaussian random variables, under such condition, 
increases with time. 


Let denote the instantaneous variance at time t^^. 


By definition * ■( exp I ( j(2C^(l+T )) 


]} = 


(A. 3.1) 


exp (2j<^) exp(-^^/2o-?) 

= J L. 


where 


PCt^) 
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The value of above integral reduces to 


2 2 

+to cos (2^) exp 12a.') 


-00 


■/2fr 


= exp (-2«yj) 


2 

For sufficiently large t the value of exp(-2<yj^) goes to zero, 
hence the value of the eq. (A. 3.1) goes to zero. 


and 
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APPENDIX - B.1 

Consider the term Re { ^ 

here 

E^c^t) = exp (jw^t + jKjP^(t) + 2jKj^P2(t) + j^j^(t)) 

and 

E2^(t) = (P2(t))^^^ exp 02t + jK2p2(t) + 2jK2Pj(t) + 
where 

P^Ct) = P^Cl + sinC^Jj^t)) with |A^ sinC^^j^Ct))! < 1 
and 

p2Ct) = p2(l + A 2 8in(to2*)) 1^2 ^ ^ I ^ 

p^ and p2 are average laser source powers, 

A^ and A2 are constants and control the depth of modulation. 

On the assumption that and are closely spaced optical 
carriers such that “ ^ 2 ’ write the term 

Re {ej^CD B;^(t)} 

. Ee|[[p^(l + Aj^ aln(eijt))]tp^(l * A^ einCUjtJJjj 

exp(- jKj^p^Aj^ ain(iOj^t) + **^*^(“2'’^^^ 

exp (j^^(t) - j<fr2Ct))| 


(B.l.l) 
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where p 


las.r .ourc. power of both Leer., 
is the average laser 

By taking approximation that 

= (l+x/2) for 1*1 < 1 

1 I'i written es 

envelope of the eq. (B.l.D o 

= CltCh,/2) elnCw.tll (ltCh,/2) 

. / + A A /SlcosC (i*!- 

- coaCCWj^tei^ltlj 


The 


(B.1.2) 


(B.1.3) 


and the phase part of eq. (B 


.1.1) ia given by [6] 


exp {-J(2isin(e.jt5 t iPjOlnCw^t) t J^^Ct) - 


=IP 

n 1 


^ (B.1.4) 


rt - K t> Aa and fly ~ *^i**o^2‘ 

where - *“ 1^0 1 ^ 

The following identities ere known. 


^1 I 

(A^/2) sinC^.>j^t) - 4j j 

’exp(j^^lt) - exp(->it)^ 

L 

i 

^2 ' 

(A^/Z) sin(402t) = XI 

■1 

[^exp(j^2^^ ■ 

1 

h 

(A 1 A 2 / 2 ) " 8 


expCKwr“2)« t expC-3«Oi-“25«] 

^1^2 

CAeA,/2) cos((«e+‘^2^*^ 8 

[expCjCe.it<Oj)t) t expC-J«0i*“2>‘>] 
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Substitutinft the above identities in eq . (B.1.1) and using the 

resut of eq. (B.1.4) in eq. (B.1.1) we get the following 
expression 


2 2 ^n^^l^ exp Cjn<^^t + jl^o^t + j^^(t) 

n 1 






n 1 


-A. 


+ P, 


• * 

Hit 


n 1 


n 1 


2 Y V ^2 

^ *^o 2 2 TT 
n 1 


Jn(/’i)Ji(/’2^^'^^*^ exp( j«i(n+l)t+jl<02*+J^l(t)--i^2^^^^ 

^n^^l exp(j<^j(n-l)t+jl^^2t+i^l(t 5-3^2^'*^^^ 

exp(>j^nt + j(l + l)6>2t+j^jCt)-j^2(^3) 


^n^^l^^l^^ 2 ^^"^^” exp( j«int+j(l-l)oi2t+j^l(t)-j^2<^3 5 


n 1 


. 2 Y Y ^1^2 

^ *^o 2 2 -T- 
n 1 


expOi(n+l)t + j(l-l)6>2t 

+j^l(t)-j^2^^>) 

JnC^l)JlC/’2^^’^^*' exp(j«i(n-l)t+j(l+l)«2t 

+j^l(t)-j^2^^^^ 


^ -I 

n 1 

n 1 


expC j<.jj (n+l)t+ j(l + l)<.> 2 t 

exp( j6>j^ (n-l)t+ j(l-l)<.> 2 '*= 


(B.I.! 
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Coll.ctlng the Ilk. frequency terms in the ebove 
eimplify it to the following expression. 




w t 




ke{pn 1 ‘"-I 

n 1 


+ jB ,) expCjnw^t + jl<*>2'*^ "*■ 




where 

^1 = ^n‘'’l 




Aj^A 

+ 8" "n+l 


J .q(/^ 


n+1 ^1^2 - frt 1 Jq 

l) 5" 1-1 2 


n 


AqA^ V, 


and 


®nl " " 4^ 


4 

A, 


_1)*^~1 - 


(i ^2 u - ^n+1 


V'’l’ •’l-l«'’2>«-^>" - r 'nf'’l> 

Equation (B.1.65 can be compactly vritt.n as the folloulng 

“‘K I 1 ‘^nl * ^^“2^ * ^2* ' 

V. *1 


n 1 


where * ^®nl 
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APPENDIX - B.2 

To compute the Fourier traaaform of 


1/2 Re 


{22 ^nl (-I’-l 2n(Ai>j+Ai^2) 4 

n 1 


(B.2.1) 


Ue know the following results 


a) Re(A^) 


1 1 * 
J ^ 7 


b) j^exp(-|T| 2iT(AVj^+Av2))j 


'[ 


2 (fiAw ) 


-1 




Aw J 


where Aw = Aw^ + Aw„ 

1 2 

c) If G(f) is Fourier transform of g(t) then 

^^l^expC ju^t) g(t)j = G(<.> - ) 

d) -T-^j^ACt) + g(t)j = ;r^j^A(t)j +^^[^fi(t)J 


Applying the result given in (a) to eq. (B.2.1) we get 

^nl ^nl ZnCAWj^+Aw^) ) 

^ n 1 

+ complex conjugate! 
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By applying the resulta given in (b), (c) and (d) in eq. (B.2.1) 
for computing the Fourier tranaform, 
ve get 
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